Draft version November 1, 2012 

Preprint typeset using I^TeX style emulateapj v. 12/16/11 



PROPER MOTIONS AND ORIGINS OF SGR 1806-20 AND SGR 1900+14 

Shriharsh p. Tendulkar 

California Institute of Technology and 
1200 E California Blvd, MC 249-17, Pasadena, CA 91125, USA 



P. Brian Cameron 

The Aerospace Corporation and 
15049 Conference Center Drive, Chantilly, VA 20151-3824, USA 



Shrinivas R. Kulkarni 

California Institute of Technology and 
1200 E California Blvd, MC 249-17, Pasadena, CA 91125, USA 
Draft version November 1, 2012 

ABSTRACT 

We present results from high-resolution infrared observations of magnetars SGR 1806— 20 and 
SGR 1900+14 over 5 years using laser-supported adaptive optics at the 10-m Keck Observatory. Our 
measurements of the proper motions of these magnetars provide robust links between magnetars and 
their progenitors and provide age estimates for magnetars. At the measured distances of their puta- 
tive associations, we measure the linear transverse velocity of SGR 1806— 20 to be 350 ± lOOkms"^ 
and of SGR 1900+14 to be 130 ±30 km s~^. The transverse velocity vectors for both magnetars point 
away from the clusters of massive stars, solidifying their proposed associations. Assuming that the 
magnetars were born in the clusters, we can estimate the braking index to be ~ 1.8 for SGR 1806— 20 
and ~ 1.2 for SGR 1900+14. This is significantly lower than the canonical value of n = 3 predicted 
by the magnetic dipole spin-down suggesting an alternative source of dissipation such as twisted 
magnetospheres or particle winds. 

Subject headings: MsignetdiTs: Neutron Stars: SGR 1806-20: SGR 1900+14 



1. INTRODUCTION 



Magnetars were proposed ( [Thompson fc Duncan|1995 



|T996) as a unified model to explain the phenomena o: 
soft gamma repeaters (SGRs) and anomalous X-ray pul- 
sars (AXPs). Magnetars, unlike canonical radio pul- 
sars, would have a very high magnetic field strength B 
(~ 10^^ G) such that their internal energy was dominated 
by their magnetic energy rather than their rotational en- 
ergy. The SGR fiares were explained as resulting from 
violent magnetic reconnections and crustal quakes and 
the quiescent X-ray emission of AXPs (which is much 
larger than their spin-down luminosity) was attributed 
to the decay of intense magnetic fields. T he discovery 
of large period derivatives (P ~ 10~^^ s s~ ^ ; [Kouveliotou' 
[et al.,,1998 ) confirmed the basic expectation of the mag- 
netar model. For recent reviews of observational and 
theoretical progress in the field we r efer the readers to 
Mereghetti] p008) and [Hurley | ( |201lD . 



Despite the successes of the magnet ar model, we have 
little understanding of why only some neutron stars are 
born as magnetars. Originally, [Thompson fc Duncan 
( [1993[ ) invoked a rapidly spinning one to three mil- 
lisecond) proto-neutron star as essential for strong am- 
plification of a seed magnetic field. The rapidly spinning 
neutron stars would result in a supernova more energetic 
than a canonical core-collapse supernova. 
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The observational support for the formatio n mech- 
anism of mag netars appears to be lacking. 



Vink 



Kuiper (2006) showed that the three supernova rem- 
nants (SJNKsj to which three magnetars are best paired, 
(Kes73, CTB 109 and N49), are completely consistent 
with the standard supernova explosion energies. 

The offset between SGR 0525— 66 (previously known as 
"5 March 1979") and its surrounding supernova remnant 
N49 and the notion that some halo SGRs might explain a 
fraction of GRBs led t o the expectation of SGRs havin g 
high space motion (see Rothschild fc Lingenfelte r|1996 ). 
This spawned a number of efforts to measure the space 
motions of magnetars. 

Here, we present astrometric observations of two of the 
youngest magnetars: SGR 1806-20 and SGR 1900+14. 
The resulting measurements of proper motion allow us to 
trace back these two objects to their potential birth sites 
and additionally measure the space motions as well. The 
paper is organized as follows. In Section [2j we summarize 
our knowledge of these two magnetars. In Section [3j we 
describe our observations, data reduction methodology 
and analysis techniques for point spread function (PSF) 
fitting, relative astrometry and photometry. We present 
the results in Sections [4] and in Section [5] we discuss the 
significance of our proper motion measurements. 



2. TARGETS 

Table □ summarizes the essential characteristics of 
both our targets; SGR 1806-20 and SGR 1900+14. We 
discuss each target in further detail in the following sec- 
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TABLE 1 

Characteristics of SGR 1806-20 and 
SGR 1900+14. 

SGR 1806-20 SGR 1900+14 



Period P (sec) 



p(io- 



P/P (kyr) 

Bsurf (1014 G) 

R.A (J2000) 
Dec (J2000) 



7.6022(7) 
49 
0.32 

24 

18^ 08™ 39.337^ 
20° 24' 39.85'' 



5.19987(7) 
17 
1.8 
7.0 

IQh 07m ;l4.3is 
9° 19' 19.74" 



Note. — Refer to 

http:/ /www. physics. mcgill.ca/^pu lsar/ ~| 
|magnetar/ main. html Positions are trom (Jtiandra 
X-ray observations. 

^ Average period derivative calculated from X-r ay p e- 
riod measurements from literature. See Section [5721 



tions. 



2.1. SGR 1806-20 



SGR 1806-20 (previously known as GB790107) was 
identified as a repeating ^ amma-ray burst witii a soft 
spectrum by Laros et al. (11986). SGR1806— 20 is best 
known for i t s giant burst of Dec ember 27, 2004 (Hurley 



et al. 2005 Palmer et al. 2005) which was one of the 



brightest cosmic hares ever detected. T he burst was fol- 



lowed by a long lived radio afterglow ([Cameron et al. 
2005[ [Gaensler et~ar]|2005[ ISpreeuw et al.||2010| ) wEicIi 
allowed the precise localization of the source . 

2.1.1. Association with Star Cluster 



SGR 1806- 20 lies in a radio nebula GlO.0-0.3 ( [Kulka- 
rni et al.|1995 ) which is a part of the W31 HIT complex. It 
was earlier suggested that the massive star LBV 1806—20 
and its surrou nding radio nebula were ass ociated with 



SGR 1806-20 



1^951) but precise 
|l)02) proved that 



( van Kerkwijk et al. 
Chandra localization ( [Kaplan et al. 
SGR 1806-20 was 14'' away trom the center of GlO.0-0.3 
and 12'' away from LBV 1806—20. A cluster of massive 
stars, c oincident with a m id-IR nebulosity, was discov- 
ered by Fuchs et al. (1999) about 7" to the north of the 



magnet ar. 

Table [2] lists all the distance measurements reported to 
date. We place a higher premium for distance estimates 
related to the X-ray counterpart of SGR 1806— 20 or the 
associated cluster of massive stars over the estimates to 
LBV 1806-20, since it is unclear whether LBV 1806-20 
is physically near the magnet ar. In Table |2j measure- 
ments 1-4 are distances to SGR 1806—20 or the cluster of 
massive stars and measurements 5 and 6 are distances to 
LBV 1806—20. We adopt a nominal distance of 9 ± 2 kpc 
which is consistent with all the measurements. 

2.1.2. IR Counterpart 

Figure [1] shows a 2 x 2 arcsec cutout near SGR 1806—20 
from our laser guide star adaptive optics (LGS-AO) sup- 
ported observations in the Kg band using the NIRC2 in- 
strument (See Section [3] for details). Star A was sug- 
gested as the NIR counterpart for SGR 180 6-20 by [Ko^ 
sugi et al.| (|2005| and independently by Israel et a 
based on NIR variability over the 2004 active 




Fig. 1 . a 2 X 2 arcsec cutout near SGR 1806-20 from a Kg band 
LGS-AO supported observation from the NIRC2 camera. Th e IR 
counterpart, as identified by (Kosugi et al. 2005 Israel et al. 2005|> 
is marked with cross hairs and labeled A as per ( |Israel et al.|2nn31 > 
as are stars B and C. 

period. Using the NAOS-C ONICA instrument on the 
8.1-m Very Large Telescope, Israel et al. (2005| moni- 
tored SGR 1806—20 on 11 epochs between March and 
October 2004. They measured a factor of two increase 
in the flux of the star A with a > 9-<j confidence. The 
IR flux increase corresponded well with X-ray flux that 
also increased by a factor of two in the 2-10 keV and 20 



100 keV bands {XMM-Newton, INTEGRAL; [Mereghetti 
et al.|[2Q05b|a[ ). Our photometric measurements show a 



factor of three varia bility in the brightness of the same 
object (Section |4.1[ ). The identification of the IR coun- 
terpart of SGR 1806—20 appears to be secure. 

2.2. SGR 1900^14 

The first bursts from SGR 190 0+14 (origina lly known 
as B1900+14) were identified by |Mazets et al.| (ri979). A 
very bright flare was detected on August 27 1998 with 
a 7-ray peak followed by a 300 - s long tail (Hurley et al. 



1999[ Kouveliotou et al. 1999). Following the burst, a 



lading radi o dl^rail et al. 1999) and X-ray source (Hurley 



et al.|p!999 ) was discovered. These observations led to a 
precise localization to within 0.1 5'^ 

2.2.1. Association with Star Cluster 

SGR 1900+14 is located near two objects from which 
it could hav e originated. A cluster of massive stars ( Vrba 
et al.|[2QQQ ), hidden behind two bright M5 super-giants, 
lies 12^^ the east of SGR 1900+14 and a 10"^ yr old, 
12' diameter SNR G042 8+00.6 lies 17' to the so uth- 
east (Mazet s et al.| | 1979| [K ouveliotou et al.||1993[ [Ya^ 
sisht et al.||l994| ). If^R 1900+14 was associated with 
the cluster of massive stars then it implies a young age 
and a space velocity close to the canonical value for pul- 
sars. However, if it is associated with the supernova 
remnant then it would have a very high proper motion. 
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TABLE 2 

Distance to SGR 1806-20 measured by various authors. 



Reference 


Distance 


Comments 




(kpc) 





Cameron 



et al. | (|2005|> 

nnffifhsSTQ 



McUlure- 
Svirski et al. (2U 
Bibby et al.| i^06E\ 



.j2005|> 

h s_^ Gaensler] ( |2005| > 



Figer et al.|(|2004|) 


Eikenberry et al. 


(2004 



6.5 - 9.8 HI absorption from Dec 2004 flare 
> 6.5 HI absorption from Dec 2004 flare 

9.4 — 18.6 X-ray scattering echos 
8.7^^ 5 Spectral classiflcation, IR photometry 
and cluster isochrones 
11.8 ± 0.5 Radial Velocity (RV) of LBV 1806 
15+^ 3 RV of LBV 1806 and surrounding nebula 

and Galactic rotation curve 
strictly > 9.5 Luminosity of cluster stars 
strictly > 5.7 Ammonia absorption to LBV 1806 



An upper limit to the proper motion (based on Chandra 
X-ray observatory imaging observations) of < lOOmilli- 
arcsecondyr"^ is nominally inconsistent with the asso- 
ciation of SGR 1900+14 with the SNR ( [Kaplan et al 
[2Q09 [ | de Luca et al 2009). 



Wachter et al. "( |2008J reported the discovery of an in- 
frared elliptical ring or shell surrounding SGR 1900+14 
which was interpreted as a dust-free cavity created by 
the giant flare of August 1998. The authors concluded 
that SGR 1900+14 is unambiguously associated with the 
afore mentioned star cluster. 

With adaptive-optics assisted Keck/NIRC2 imaging 
and Keck/ NIRSPEC sp e ctrosc opy of the cluster near 
SGR 1900, Id avies et al. rt2009k estimated the progeni- 
tor mass to be 17 + 2 IVl© which is much lower than the 
progenitor masses estimated for other magnetars (~ 40 
to 50 M0). 

2.2.2. Distance 

'Vrba et al.' ("19 961) showed that the bright IR sources 
noted by Hart mann et al | ^996^ at the ROSAT lo- 



calization of SGR 19UU+14 were M5 super-giant stars 
at a distance of 12 to 15 kpc with an extinction of 
Ay ^ 19.2 mag. Davies et al.| (2009) measured a ra- 



dial velocity of —15.5 + 4kms~^ for the cluster of stars 
implying a distance of 12.5 + 1.7 kpc using the measured 
m odel of Galac t ic rot ation. We adopt the measurement 
oflDavies et al.|(|2009|) for the distance to SGR 1900+14. 



2.2.3. IR counterpart 

Figure [2] shows a 4 x 4 arcsec Kg band image from our 
LGS-AO observations with the NIRC2 camera around 
the X-ray pos ition of SGR 1900+ 14. The stars are la- 
belled as per Testa et al. (|2008|). They obtained two 
Ks band AO observations of the same field around 
SGR 1900+14 with VLT NACO instrument in March 
and July 2006. Star 7 was the only source inside the 
radio-position error circle (dashed circle in Figure^ that 
showed a photometric variability. They detectea a 3-cr 
increase in the fiux of star 7 and proposed it as the IR 
counterpart of SGR 1900+ 14. We accept the counterpart 
proposed by [Testa et°aD (|2008). 

In an attempt to gather additional evidence for the 
identification of the IR counterpart, we have measured 
Kp band photometric variability and H — Kp colour for 
the stars in the field. These measurements are reported 




Fig. 2. — A 4 X 4 arcsec cutout near SGR 1900+14 from a Kg band 
LGS-AO supported observation from the NIRC2 camera. Stars are 
labelled as per Testa et al. (2008). The black circle is centered on 
the radio position of yGK19UU+14 from Frail et a l. (1999) and 
encircles the 0.8^^ radius, 99%-confidence circle trom Testa et al.| 
( 2008 ) the positions from which are used for our absolute astrom- 
etry. Star 7 is the proposed counterpart of SGR 1900+14 based on 
its variability. 



in Section |4.2[ However, during this period, the X-ray 
counterpart did not show significant variablity. Hence 
the absence of NIR variation of the proposed counterpart 
does not provide any new insights. 

We report (in Section 4| that the proper motion of 
star 6 lies along the Galactic rotation curve, whereas 
the proper motion of star 7 is significantly different 
from those of galactic stars. This evidence strength- 
ens the identification of star 7 as the IR counterpart of 
SGR 1900+14. 

3. OBSERVATIONS AND ANALYSIS 

3.1. Observations 

Starting in 2005 to the present time, we undertook a 
program for astrometric monitoring of magnetars with 
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TABLE 3 
Summary of observations of 
SGR 1806-20. 



Date & MJD 


Filt 


Cam 


Exp 


(UTC MJD) 






(s) 


2005-03-04 53433.641 


Kp 


N 


1440 / 


2005-04-30 53490.511 


Kp 


N 


750 


2005-08-10 53592.366 


Kp 


W 


600 


2005-08-11 53593.344 


Kn 


w 


840 


2005-09-26 53639.258 


Kp 


N 


600 / 


2006-07-03 53919.403 


Kp 


N 


2820 


2006-08-17 53964.304 


Kp 


N 


1800 / 


2007-05-22 54242.487 


Kp 


N 


1020 


2007-06-11 54262.403 


Kp 


N 


2040 


2007-07-16 54297.345 


Kp 


N 


3000 


2007-08-06 54318.329 


Kp 


N 


2640 / 


2008-05-21 54607.468 


Kp 


N 


2460 


2008-06-29 54646.407 


Kp 


N 


3360 


2008-07-26 54673.342 


Kp 


N 


3180 / 


2010-06-18 55365.442 


Kp 


W 


80 



Note. — A / in Column 4 marks the im- 
ages used for astrometric measurements. 



TABLE 4 
Summary of observations of 
SGR 1900+14. 



Date & Time 


Filt 


Cam 


Exp 


(UTC) 






(s) 


2005-04-30 53490.558 


Kp 


N 


1300 / 


2005-08-09 53591.434 


Kp 


W 


2400 


2005-08-10 53592.400 


Kp 


W 


300 


2005-09-26 53639.349 


Kp 


w 


720 


2006-07-03 53919.472 


Kp 


N 


1980 


2006-07-04 53920.511 


Kp 


N 


1920 


2006-08-17 53964.439 


Kp 


N 


1140 / 


2006-10-13 54021.242 


Kp 


N 


2220 


2007-05-22 54242.550 


Kp 


N 


1500 / 


2007-06-11 54262.553 


Kp 


N 


1260 / 


2007-06-11 54262.582 


H 


N 


660 


2007-08-06 54318.455 


Kp 


N 


1800 / 


2007-11-03 54407.229 


Kp 


N 


1260 / 


2008-05-21 54607.564 


Kp 


N 


1260 


2008-06-29 54646.471 


Kp 


N 


3660 


2008-07-26 54673.405 


Kp 


N 


2280 / 


2008-10-22 54761.227 


Kp 


N 


1680 / 


2009-04-06 54927.599 


Kp 


N 


1620 


2009-07-17 55029.340 


Kp 


N 


2100 / 


2009-08-04 55047.346 


Kp 


N 


2100 


2009-09-29 55103.226 


Kp 


N 


2340 


2010-06-18 55365.470 


Kp 


N 


2340 / 



Note. — A / in Column 4 marks the im- 
ages used for astrometric measurements. 

the 10-meter Keck 2 telescope using the Laser (Suide 
Star Adaptive O ptics (LGS-AO; [Wizinowich et al.||2QQ6 
van Dam et al. 2006) and the N ear-infrared Camera 2 
(iNlRC2). The log of 
Tables El and H 



3g of our observations can be found in 



3.1.1. NIRC2 

The NIRC2 instrument has two modes: wide (W) 
and narrow (N) with a field-of-view (FoV) of ~ 10 x 
lOarcsecond and ^ 40 x 40arcsecond respectively. The 
corresponding pixel scales are 9.942 milli-arcsecond per 
pixel and 39.768 milli-arcsecond per pixel. The wide field 
images were obtained to aid transferring the photometry 



and astrometry from the low resolution 2MASS images 
to the small FoV narrow camera NIRC2 images. The 
narrow field images were used for the astrometric mea- 
surements. Based on weather and faintness of each mag- 
netar, multiple short (~20 s) exposures were chosen to 
avoid saturating the detector. The typical full width at 
half maximum (FWHM) achieved in these observations 
was ~ 70 milli-arcsecond ~ 7 pix. 

Each of the NIRC2 narrow camera images was in- 
spected for quality control. Images in which the AO 
correction was poor were rejected. The shallow images 
with acceptable AO correction were rejected for astrom- 
etry due to the non-detection of the magnetar and/or 
lack of sufficient reference stars but were used to pho- 
tometric calculate upper limits on the brightness. The 
images used in the final proper motion measurement are 
denoted by a / in Column 4 of Tables |3] and [4j 

3.2. Data Analysis 

The images from the NIRC2 camera were reduced 
using the FITS analysis package pyraf in a standard 
manner by subtracting corresponding dark frames and 
flat-fielded using appropriate dome-flats. A sky fringe 
frame was made by combining dithered images of mul- 
tiple targets with the bright stars m asked. We used 
SExtractor (Bertin & Arnouts 1996) for the prelimi- 
nary detection and masking of stars. The fringe frame 
was subtracted after being scaled to the appropriate 
sky background level. Before coadding the frames, each 
frame was corrected for optical distortion using a distor- 
tion solution measured for NIRC^ 

3.2.1. PSF Fitting 



We used the IDL package StarFinder ( Diolaiti et al. 
2000) to perform PSF estimation, fltting and subtrac- 



tion. This code iteratively estimates a normalized PSF 
shape from user selected stars, while subtracting faint 
neighboring stars to minimize the contamination of the 
PSF estimate. StarFinder flts a constant PSF shape 
over the entire fleld of view (FoV). This assumption ap- 
pears to work well for the NIRC2 narrow camera FoV. 
The uniformity of the PSF over the FoV also mitigates 
the errors from centroiding variable PSFs. 

AO PSFs differ from PSFs obtained from atmospheric 
seeing limited observations in two aspects: Firstly, be- 
cause the AO correction decorrelates as a function of 
distance from the AO reference source (i.e. sodium laser 
beacon), the PSF varies radially across the field of view. 
Secondly, since AO correction cannot correct all of the 
wavefront errors caused by atmospheric turbulence, even 
on- axis, AO PSFs have a distinctive shape with a sharp 
diffraction-limited (FWHM X/Dtei) core and a wide 
(FWHM ~ atmospheric seeing) shallow halo around it. 
For the Keck AO system, these components are 44 milli- 
arcsecond and ~ 1 arcsecond respectively. The order of 
magnitude difference in size and brightness of the two 
components makes it challenging to accurately measure 
and subtract the PSF in the image. We describe how 
both these challenges are handled in the next paragraph. 

To further reduce the effect of PSF variations, relative 
photometry and astrometry measurements were down- 



^ See http://www2.keck.hawaii.edu/inst/nirc2/ 
forReDoc /post _o bser ving / dewarp / ^ 
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weighted farther away from the object under considera- 
tion. The detail s of the relative weighting are described 
in Section 13.2.21 The PSF model size was chosen to be 
200 pixels (1.95 arcsecond) wide to encompass both the 
core and the halo of the PSF. The few brightest stars 
in each of the fields were used for estimating the halo 
contribution. 



3.2.2. Relative Astrometry 

demonstrated a framework for 



Cameron et al. 



(|2009| 

high precision astrometry (< 100 /iarcsecond) through 
an optimal estimation technique that availed the correla- 
tions in stellar position jitter. We use the same method- 
ology with modifications for including the proper mo- 
tions of the stars over multiple epochs and an appropriate 
weighting scheme. 

The dominant source of astrom etric error in the single 



epoch, short exposure images of Cameron et al. (2009) 



was tip-tilt anisoplanatism. For our coadded long expo 
sure images the tip-tilt anisoplanatism is averaged out. 
We constructed the covariance matrix theoretically us- 
ing geometry of the field and a typical turbulence pro- 
file from Mauna Kea. The residual distortion of the 
NIRC2 distortion solution has a root-mean-square value 
of 1 milli-arcsecond. However the distortion residuals 
have higher values towards the edgeaj To reduce the 
effect of residual distortion, especially in images with 
significant dithering, a separation- weighted measurement 
scheme (the term used below) was used to downweight 
stars far from the target. 

To account for the proper motions of all the stars in 
the field, it was necessary to include the proper motion 
estimates in the framework and simultaneously estimate 
a least-squares fit for grid positions and proper motions. 
Given N stars detected in the field, the measurement 
of the offset between the target star and each of the re- 
maining stars results in a set of vectors at each of the E 
epochs. 

The differential offsets between star and the grid of 
N reference stars at epoch k is written as a single column 
vector, 

do/c = [^01, • • • ,^OAr,^oi, • • • ,yoN]k' 

Here the distance between the x- 

coordinate of the j*^ reference star and the x-coordinate 
of the i^^ target star, and likewise for y. The goal of dif- 
ferential astrometry is to use d to determine the position 
of the target star with respect to the reference grid of 
stars at each epoch. 

We use a linear combination of the elements of d with 
weights Wi to obtain the relative position of target star 
i at epoch /c, 

Pi/c Widi/e, 

where, for example, the weight matrix for star 0, Wq is 



Wo 



w 



yx,01 



We calculated weights as follows: 



'^xy.QN 
^yy,ON 



afj. Here afj 



'^^Tj^ij^ where <j^j is the geometric 



^ http:/ /www2. keck. hawaii.edu/inst/nirc2/ 
forReDoc / post _observing/dewarp / 



mean of the parallel and perpendicular comp onents of 
the tip-tilt j itter as defined in Equation 1 of | Cameron 



et al. (20091); and Oij is the angular offset between the 
star I and tne star j. We have used the notation Wxy,oj 
to denote the weighting of the offset from the target star 
{i = 0) to star indexed j in the y direction which is used 
to determine the x component of the target's position, 
p. 

We assume a simple linear model for the stellar motion 
where x = Zx -\- Vxt- The differential offsets are thus a 
column vector, 

' Zx,l + Vx,lt - {Zx,0 + Vx,ot) ■ 
Zx,N + VxM - {Zx,0 + Vx^ot) 

Zy^N + "^y.Nt — (%,0 + ^y,0t)- 



do 



and the unknown quantities are. 



\Zx,0^ • • • 7 Zx,N^ Vxn, . 



(1) 

(2) 



We solve for the variables b from the vector d given 
weights W in the least squares sense. For a given target, 
we use the same weights for all epochs. The overall x 
and y shifts of each image (i.e. the registration of the 
image) are fit as free parameters in this method. 

NIRC2 is mounted at the Nasmyth focus of the Keck H 
telescope. A field-rotator allows the observer to set the 
position angle of the instrument. Our default position 
angle was zero degrees (North is up and East is to the 
left on the detector). However, there are small errors in 
the setting of the field rotator as well as tracking errors. 

To measure this, we chose the images obtained on May 
22, 2007 as the reference image for both the targets. The 
reference images were chosen on the basis of good AO 
correction and image depth. We computed the rotation- 
angle and the plate-scale of the image at each epoch with 
respect to the reference image. We find that the rotation 
angle is within 0.5 degrees and the image scaling is within 
0.1% relative to those of the reference image. The stellar 
position grids were corrected for the measured rotation 
and plate-scale changes before measuring their proper 
motions. 

To understand the systematic effects caused by our 
choice of grid stars, we re-analyzed the centroiding data 
after randomly eliminating a selected number of stars 
from the reference grid. We compared the results to those 
obtained from our entire grid of stars. For example, by 
eliminating one randomly chosen star out of the 50 stars 
in the SGR 1900+14 field, the proper motions of all other 
stars change by A{fia,l^6) = (7.6 ± 15.4, 17.1 ± 13.7) x 
10~^ milli-arcsecondsyr"^. This is much smaller than 
our statistical errors of ~ 1 milli-arcsecond yr~^. Hence 
we conclude that the choice of our reference grid is robust 
and does not add significant errors to our measurements. 

3.2.3. Galactic Rotation 

Since our relative astrometry framework calculates the 
proper motion of each object with respect to a grid of 
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neighboring stars (i.e. with respect to the average mo- 
tion of ah other stars), it imphcitly assumes that the net 
velocity of the field is zero. However, this is not true 
since the rotation of the Galaxy and the peculiar veloc- 
ity of the Sun with respect to the local standard of rest 
(LSR) cause significant motions at the precision we seek. 
Our framework cannot measure the net velocity of the 
field without prior knowledge of the absolute motion of a 
few stars or equivalent ly, the absolute non-motion of an 
extra-galactic object in the field. 

To correct for this effect, we need to calculate the mean 
galactic proper motion of all the stars in the field along 
the line of sight given by Galactic longitude and latitude 
(/, b). We modelled the differential rotation of the Galaxy 
and the local velocity of the Sun and calculated the effec- 
tive proper motion of an object at a given position (r, b) 
in the Milky Way, where r is the distance away from the 
Sun. We made a model assuming the local vel ocity of the 
Sun to be (U .V.W) = (10.0, 5.2, 7.2) kms-^ (|Dehnen 



Binney|1998 ) and that the Galaxy is rotating with a con- 



stant circular speed outside of Ri = 2 kpc of 220 km s 



decreasing linearly inside of that Ri (Binney & Tremaine 
[2008) . We set the distance fro m the Sun to the cente r 
ot the Galaxy to Rq = 8.0 kpc ( Eisenhauer et al.| |2003). 



From the rotation curve, we calculate the Galactic proper 
motion /iQai = [Ma ^M^] Gal of objccts at various distances 
(Ikpc < r < 20 kpc) in the direction (/, b) of the magne- 
tar that are moving with the Galactic flow. 
We estimate the number density o f stars in the Milky 



Way using the model calculated by Juric et al. ( [2008 ) 
using SDSS data. They fit a thin disk, thick disk and 
a halo to the SDSS data set and calculate the number 
density function based on their fit. Along the line of 
sight, the number of stars in our field at a distance r from 
the Sun is proportional to r'^p{R,Z)^ where p{R^Z) is 
the number density of stars at the cylindrical coordinates 
(i?(r, /, 6), Z(r, /, b)) in the Milky Way. 

For a given field, we calculate the velocity of the field 
MFieid = [Mq;7 M(5] Field cis the integral of the proper motion 
weighted with the number density as described above. 
This gives. 



MField 



j"^^^^ r'^p(^R{r,l,b),Z{r,l,b)) x [/ic.,5(r, 6)Gai]dr 



/;r: rM^(r,/,6),Z(r,/,6))dr 



Thus, the total proper motion of each object in the 
sky is /isky,i = MR,i + /^Fieid- Table [sllists the calculated 
proper motion for the field and the ualactic proper mo- 
tion for an object at the distance of the magnet ar for 
both of the targets. 

3.2.4. Peculiar Motion 

We are interested in back-tracing the proper motion of 
the magnetar to identify its birthsite and estimate the 
time since it left the birthsite. The relevant motion for 
this measurement is the relative proper motion between 
the magnetar and its progenitor. A reasonable assump- 
tion is that the progenitor, likely a young massive star, 
was moving with the Galactic rotation curve. We define 
the peculiar motion of the magnetar as the difference 
between its total proper motion /2sky,i and its expected 
Galactic proper motion /icai, i-e. /isky,i = A^Gai + /^Pec- 



With this definition, the transverse velocity of the mag- 
netar relative to its neighbourhood becomes r|/ipec| in a 
direction ^, s.t. tan(6>) = {fJ^a/ fJ^5)pec East of North. 

3.2.5. Photometry 

StarFinder calculates flux estimates for stars in the 
field by scaling the normalized PSF model to best fit the 
image. We calculate the photometric zero-point (ZP) 
for each image by comparing the magnitudes of stars 
to the 2 Micro n All Sky S urvey (2MASS) Point Source 



Catalog (Skrutskie et al. 



2006) and to published high- 



resolution studies of the iields which were anchored to 
the 2MASS catalog. The details of comparison stars for 
each field are given in Section [4j 

4. RESULTS 

4.1. SGR 1806-20 

We performed PSF fitting on the NIRC2 narrow cam- 
era images to identify 71 stars through 10 epochs. The 
positions of these 71 stars were used for relative astrom- 
etry. 




Epochs (YY/MM) 

Fig. 3. — Kp magnitudes of stars around SGR 1806—20 measured 
over period of 3 years. The circles (red in the onhne version) corre- 
spond to star B and squares (blue in the online version) correspond 
to star C. The counterpart (star A) of SGR 1806— 20 is marked by 
black triangles. We note a clear variation over a factor of 3 in the 
brightness of star A. 

We performed relative photometry on the stars A, B 
and C in Figure[l] The photometric zeropoints were mea- 
sured by matching t he magnitudes of sta rs B and C to the 
values measured by [Kosugi et al. (2005). Figure [s] shows 
the measured magnitudes of the three stars. We observe 
a clear factor of 3 variation in the brightness of the IR 
counterpart of SGR 1806— 20, star A, thus securing the 
identification of the IR counterpart of SGR 1806— 20. 

4.1.1. Proper Motion 

Figure |4] shows the measured proper motions of the 
stars in the SGR 1806—20 field. The field velocity correc- 
tion was calculated to be (/ia, /i^) Field = (3.0, 4.8) milli- 
arcsecondyr"^. The proper motion of SGR 1806— 20 
away from a putative progenitor in the galactic fiow is 
(/icK, /i^) = (—4.5 ± 1.4, —6.9 ± 2.0) milli-arcsecondyr"^. 
Assuming a distance of 9 ± 2 kpc, this corresponds to 
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TABLE 5 

Proper motions calculated from the G alacti c rotation model as described in 

Section [3^3] 



Object ID 


Distance 




MField 












[Ma,M(5] 




(kpc) 


(deg) 


(milli-arcsecond yr~ -"^ ) 


(milli-arcsecond yr~ -"^ ) 


SGR 1806-20 


9±2 


(10.0,-0.2) 


[3.0,4.8] 


[4.2 + 0.9,7.0 + 1.8] 


SGR 1900+14 


12.5 ± 1.7 


(43.0, +0.8) 


[2.7,4.6] 


[2.7 + 0.2,4.8 + 0.4] 




2 4 

/i„(mas/yr) 

Fig. 4. — The proper motion of 71 stars in the field of 
SGR 1806-20 in the sky frame of reference. SGR 1806-20 is 
marked by the star with error bars (colored red in the online 
version). The remaining stars have only their best-fit values 
(hollow black circles) after adding the bulk motion of the field 
(/^Field = (3.0, 4.8) milli-arcsecond yr"-*^) (marked by a black '+'). 
The thick gray line represents the expected motion of stars from 
1 to 22.8 kpc along this lin e of si ght, as per the Galactic rotation 
model presented in Section [3.2.3] Black dashes along the line de- 
note positions 1, 5, 10, 15 and 20 kpc away from the Sun. The 
section of the line representing objects at a distance of 9 + 2 kpc 
from the Sun is marked with a black star and black line to de- 
note the possible motion of the progenitor of SGR 1806— 20. The 
dashed diagonal line (green in the online version) is the locus of 
objects with /i^ = 0, i.e. with zero proper motion along galactic 
latitude. Other high proper-motion objects, probably halo stars 
are marked by diamonds. The square marks the nominally high 
proper-motion object near the edge of the detector. However, this 
measurement may be corrupted by distortion residuals and hence 
is not considered any further. 

a linear velocity of 350 ± lOOkms"^ with an angle of 
213° ± 10° East of North. 

Figure [5] shows the direction of motion of SGR 1806—20 
with respect to its neighbors. Backtracing this space 
velocity would put the magnetar close to the cluster of 
massive stars about 650 years ago. 

4.1.2. Other High Proper-Motion Stars 

In Figure ffl we mark the high proper- mot ion objects 
with diamonds and squares. These stars deviate signif- 
icantly from the dashed green line marking the locus of 
objects with /i^ = 0, i.e. with zero proper motion along 
the galactic latitude. These are probably halo stars mov- 
ing at a high speed through the Galactic disk. 




Fig. 5. The position of SGR 1806-20 (diamond, blue in the 
online version) traced back by 0.65 kyr is marked by the ellipse 
(colored red in the online version). The size of the ellipse de- 
notes the positional uncertainty corresponding to the uncertainty 
in the proper motion measurement. The solid lines (red in the on- 
line version) represent the 1-a limits on the angle of motion. The 
dashed circle (cyan in the online version) denotes the clust er of 
massive stars corresponding to the mid-IR source of Fuchs e t al.| 
(1999) . The position of the luminous blue variable LBV 1806^^ 
is marked. 



4.2. SGR 1900^14 

We observed SGR 1900+14 at 13 epochs with an ex- 
posure time of about 1 hour at each observation. Using 
Kp-band photometry and H — Kp band color ( at a single 
epoch), we present variability and color measurements 
of SGR 1900+14 and its surrounding stars. Our ab solute 
astrometry is matched to positions as reported by [ Testa 
|et al.| (|2008) with an accuracy of 6 milli-arcsecond. They 
reported a 3-cr astrometric uncertainty of O.SV which we 
adopt for comparison with the radio position for Figure [2] 
In three images of the SGR 1900+14 field that had ex- 
cellent AO correction, we detected a faint source (labelled 
10 in FiCTre [2| 0.2^^ away froni star 3. Source 10 is not 

([^08') as it was blended with 



detected by ITesta et al. 



star 3. However, we detected no variation in the com- 
bined brightness of star 3 and 10 in our data and the 
measurements from Testa et al. (2008) within 0.07 mag. 
Star 10 is a factor 01 ~ 40 fainter than star 3. With this 
ratio, assuming no variation in the light from star 3, we 
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TABLE 6 

Persistent X-ray luminosity of S GR 1900+14 in 

THE 1-10 KEV band as REPORTED BY IMeREGHETtFI 

I ET al7| ( [2006J AND [Israel et al.| ( |2008| )" 



Interval 


Fx 


(UTC Date) 


(10-12 erg cm-^s-i) 



20 Sep 2005 - 


22 Sep 2005 


4.8 ±0.2^ 


25 Mar 2006 - 


27 Mar 2006 


4.6 ±0.8^ 


28 Mar 2006 - 


28 Mar 2006 


6.3 ± 1.7^ 


01 Apr 2006 - 


01 Apr 2006 


5.5 ±0.4^ 


08 Apr 2006 - 


10 Apr 2006 


5.0 ± 1.4^ 


11 Apr 2006 - 


15 Apr 2006 


5.0 ±0.7^ 



^ Absorbed 0.8-12 keV flux from [Mereghetti et aT] 
(2006). 



Unabsorbed 1-10 keV flux from 



Israel et al. 



(2008]). 



can constrain the maximum variation in the brightness 
of star 10 to be 0.4 mag as compared to the 0.48 mag 
variation measured for star 7 and no variation for star 



3 reported by Testa et al. (2008). Thus, we continue to 
accept star 7 as the IK counterpart of SGR 1900+14. 

4.2.1. Variability 



17.0 
17.5 
18.0 
18.5 
' 19.0 
19.5 
20.0 
20.5 
21.0 



in ID ID ID 



00 00 00 ai 



Epochs (YYYY/MM) 

Fig. 6. — Relative photometry light-curves of stars 2-7 (except 
5) around SGR 1900+14. To reduce the effect of PSF variations 
over the field, relative photometry was performed on nearby stars 
and the absolute calibration was performed by matching stars 2, 3 
and 4 to their magnitudes as measured by Testa et al. (2008 ). The 
inverted triangles mark 3-a upper limits for star Y when it was not 
detected at the edge of star 3. 

Figure |6] shows the photometry of stars 2-7 (except 
5^^ The median magnitude offsets of stars 2, 3 and 4 
were used as relative ZP offsets and the absolute ZP off- 
sets were ca l culate d using Kp magnitudes as repor ted by 



Testa et al 
et al.|(:j008 



(2008). The counterpart suggested by Testa 

, star 7, was not detected at the edge of star 

3 on epochs when the images were not sufficiently deep 
or the AO performance was not satisfactory. The non- 
detections were marked with the upper limit on the flux 

^ Star 5 is excepted from all further discussion since it is far 
away from the X-ray position error circle and does not affect any 
of the conclusions. Its identification in the middle of the numbering 
range is an unfortunate quirk of the numbering scheme that was 
implemented in previous literature. 



(black triangles). Including the upper limits on flux, star 
7 shows slight variability but it is not conclusive. 

During our entire observation period from 2005 to 
2010, the X-ray counterpart of SGR 1900+14 showed 
burst activity in only one period from March to June 
2006 (Israel et al. 2008). Unfortunately, we have no IR 
observations between September 2005 and July 2006. Of 
these, the AO performance in July 2006 was not satis- 
factory leading to poor photometry and source confusion. 
As shown in Table |6l the persistent X-ray luminosit y as 
measured by IsraeFet al. (2008) and Mereghe tti et al.| 
([2006) showed a slight increase in March 2006 and de- 
creased to the pre-burst value by April 2006. Thus the 
lack of NIR variability is not surprising. 

4.2.2. Color Measurement 

During the June 11, 2007 observations, we obtained 
Kp and H band images of the field. These images were 
used to determine the colors of stars near SGR 1900+14. 
No high-resolution H band photometry of this field has 
been performed previously, so we chose to use 2MASS 
measurements of bright stars to calculate the ZP offsets 
for the H band image. The problem with this imple- 
mentation was that stars bright enough to be included 
in the 2MASS catalog were were saturated in the NIRC2 
images which were intended to image the faint magne- 
tar. We rely on the reconstruction of the saturated cores 
of bright stars by StarFinder. This increases the error 
in photometric measurement and hence in the ZP esti- 
mate. We estimate this systematic error in H band ZP 
to be 0.5 mag. This systematic error changes the scalin; 
on the X-axis of the color-magnitude diagram (Figure 
and should not change the conclusion if the magnetax 
were to have a color distinctly different from other stars 
in the field. 




Fig. 7. — H — Kp color vs Kp magnitude diagram for 50 stars 
in the SGR 1900+14 field. Stars 2-7 (except 5) are marked.The H 
band image zero-point has a systematic uncertainty of ~ 0.5 mag 
which would effectively only change the scale of the x— axis. 



Figure W\ shows an H — Kp color vs. Kp magnitude 
diagram for the 50 stars in the field. Stars 2-7 are la- 
beled. Neither star 6 nor star 7 have abnormal colors 
and neither is distinctive. There is no clear structure 
(for example, a main-sequence) in the color- magnitude 
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TABLE 7 

H AND Kp BAND PHOTOMETRY FOR STARS 2 7 (EXCEPT 5) NEAR 

SGR 1900+14. The zero-point error in the photometry is 0.5 mag for 

BAND. 



Object ID RA (J2000) Dec (J2000) H band Kp band 

(deg) (deg) (mag) (mag) 



2 19^07^14.28^ 9° 19' 18.84'' 18.57 ±0.003 17.98 ± 0.002 

3 19^07^14.30^ 9° 19' 19.63" 17.76 ± 0.002 17.19 ±0.001 

4 19^07^14.28^ 9° 19' 19.78" 18.96 ± 0.005 18.41 ± 0.003 

6 19^07^14.34^ 9° 19' 19.92" 20.41 ± 0.02 19.74 ± 0.01 

7 19^07^14.31^ 9° 19' 19.74" 21.17 ±0.04 20.63 ± 0.02 



diagram. This is probably due to the varied distances, 
ages and extinctions to the stars in this direction. Table[7| 
hsts the H and Kp band magnitudes of stars 2-7 (except 
5) as shown in Figure ItJ Magnetars are not known to fall 
in a specific color bana and our lack of understanding of 
the background physics prevents us from p redicting the 
shape of the IR emission spectrum (jTesta et al. 2008). 
We conclude that the lack of a distinctive color for any 
star near the location of SGR 1900+14 is not significant. 



Figure |8] shows the measured proper motions of 50 
stars in the neighborhood of SGR 1900+14. The ve- 
locity offset, calculated from the galactic rotation, is 
(/io;, /i^) Field = (2.7, 4.6) milli-arcsecondyr"-*^. For star 
7, we calculate a proper motion of {jia-, l-^d) = (—2.1 ± 
0.4, 0.6 ± 0.5) milli-arcsecondyr"^ away from a putative 
progenitor moving with the galactic fiow. At a distance 
of 12.5 ± 1.7 kpc, this corresponds to a transverse space 
velocity of 130 ± 30kms~^. 



4.2.3. Proper Motion 



^ 4 



-1 



+ 




1 kpc 



2 3 4 

/i„(mas/yr) 



Fig. 8. — The proper motion of 50 stars in the field of 
SGR 1900+14 in the sky frame of reference. The putative counter- 
part of SGR 1900+14 is marked by the star with error bars (colored 
red in the online version). The proper motions of star 6 (solid black 
triangle) and star 3 (inverted black triangle) seem to lie along the 
Galactic rotation curve. The remaining stars have only their best- 
fit values (hollow black circles) after adding the bulk motion of the 
field (/iField = (2.7, 4.6) milli-arcsecondyr"-*^) (marked by a black 
+). The thick gray line represents the expected motion of stars 
from 1 to 19.8 kpc along this lin e of sig ht, as per the Galactic rota- 
tion model presented in Section [3.2.31 Black dashes along the line 
denote positions 1, 5, 10, 15 and 20 kpc away from the Sun. The 
section of the line representing objects at a distance of 12.5+1.7 kpc 
from the Sun is marked with a black star and a black line to denote 
the possible motion of the progenitor of SGR 1900+14. The dashed 
diagonal line (green in the online version) is the locus of objects 
with jiij = 0, i.e. with zero proper motion along galactic latitude. 




Fig. 9. — The position of the putative counterpart of 
SGR 1900+14 (blue diamond) traced back by 6kyr is marked by 
the solid ellipse (red in the online version). The size of the ellipse 
denotes the positional uncertainty corresponding to the uncertainty 
in the proper motion measurement. The solid (red) lines represent 
the 1-cr limits on the angle of motion. The dashed circle fcva n in 
the online version) denotes the cluster of massive stars ( jVrba et al.| 
[2000t . 

Figure [9] shows the direction of motion of SGR 1900+14 
with respect to its neighbors. Backtracing this space 
velocity would put the magnetar close to the cluster of 
massive stars about 6kyr ago. 

Star 6 and star 3 are the only two other sources de- 
tected inside the 3-<j error circle around the radio po- 
sition of SGR 1900+14. Their velocities are marked by 
a black triangle (Star 6) and an inverted black trian- 
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TABLE 8 

Proper motions measured for stars 2 7 near 
SGR 1900+14. The values have been corrected for 

THE GALACTIC ROTATION OFFSETS. ThE TRANSVERSE 
SPACE VELOCITIES ARE CALCULATED ASSUMING A DISTANCE 
OF 12.5 kpc. 1-a ERROR BARS ON /ipec ARE 
(0.4, 0.5) MILLI-ARCSECOND YR~^ . 



Object 


PPec 


Velocity 


Direction 




(milli-arcsecond yr~ -"^ ) 


(kms-i) 


E of N 


2 


(-0.11,-0.55) 


33 + 25 


191 + 143 


3 


(-0.08,-0.67) 


40 + 25 




4 


(-0.74,-2.39) 


148 + 30 


197 + 10 


6 


(+0.88, +1.58) 


107 + 30 


30 + 12 


7 


(-2.11,-0.61) 


130 + 30 


254 + 10 



gle (Star 3) in Figure [sj Their velocities suggest that 
these are regular galactic stars moving in the plane of 
the galaxy (dashed green line). 

Table |8] gives the proper motions measured for each of 
the stars 2-7 along with their corresponding transverse 
space velocity assuming a distance of 12.5 kpc. 

5. DISCUSSION 

Using LGS adaptive-optics supported near-IR 
observations, we have measured the proper mo- 
tions of SGR 1806-20 and SGR 1900+14 to be 
{jJiai IJ^s) = (— 4.5, — 6.9) ± (1.4, 2.0) milli-arcsecondyr"^ 
and (/ia,M5) = (-2.1,-0.6) ± (0.4, 0.5) milli- 
arcsecond yr~^ respectively. These correspond to a 
linear transverse velocity of 350 ± lOOkms"^ and 
130 ± 30kms~^ respectively at the measured distances 
of their putative associations. Previously, using Very 
Long Baseline Interferometry (VLBI) at radio wave- 
lengths, transverse linear velocities have been measured 
only for two mag netars: the AXP IE 1810-197: 
212 ± 35 kms 
AXPPSRJ1550 



dHelfand et al.| |2007|) and the 



(Deller et al 



5418: 280 ±120 kms 
20T2I. The radio counterpart fo r AXP PSR' J1622- 495iy 
has been recently identified by Levin et al. ([2010) and 



would lead to an accurate proper motion measurement 
with VLBI. With the transverse velocity measurements 
for two AXPs and two SGRs in the 100 - 400 km s"^ 
range, it is highly unlikely that each of these objects has 
an extremely high radial velocity component. Hence we 
conclude that magnetars as a family do not possess the 
high space velocities (~ 1000 km s~^) that were expected 
earlier (cf. Rothschild & Lingenfelter 1996). 

Consider the space velocities of other families of neu- 
tron stars in contrast with magnetars. Canonical ra- 
dio pulsars {B ~ 10^ ^ G) have typical space velocities 
of - 200 -300 km s"^ { jHobbs et al.|2005[ ). [Tetzlaff et aL 
( [20 10) traced the motions of 4 young, hot X-ray bright 
isolated neutron stars to associate them with progenitors 
and constrain their ages. They calculated the space ve- 
locities of these objects to be 350 ± 180kms~^. There 
are a few fast moving pulsars such as PSR J1357— 6429, 
which is a Vela-like radio pulsar has a transverse velocity 
of 1600 -2000 kms- ^ dKirichenko et aH2012| , but these 
seem to be outliers from the family. From these data, 
we observe that perhaps velocities are not a good dis- 
criminator of different groups of neutron stars and their 
origins. 



5.1. Association 

Our measured proper motions provide very good ev- 
idence linking SGR 1806— 20 to the cluster of massive 
stars. The time required for SGR 1806— 20 to move from 
the cluster to its current position is 650 ± 300 yr. It 
may not be a surprise that one of the younger super- 
novae in our galaxy resulted from the magnetar. How- 
ever, SGR 1806— 20 lies in the galactic plane behind dust 
clouds which create very high extinction in the visible 
wavelengths. Hence, the supernova associated with the 
magnetar may not have been visible to the naked eye. For 
SGR 1900+14, we rule out any association with the su- 
pernova remnant G 42.8+0.6 and confirm that this mag- 
netar is associated with the star cluster. The time to 
trace the magnetar back to the cluster is 6 ± 1.8 kyr. 

The turn-off masses for the clusters with which 
the magnetars are associated allow us to place lower 
limits on the progenitor masses of these magnetars. 
Currently, progenitor mass estimates exist for three of 

the magnetars: 

SGR 1806-20: 48tf ( |Bibby et al.|2008|) , 
CXO J1647 -455: > 401 " " 

et al.||2 010') and 

SUR 1900+14: 17±2M0 ( |Davies et al.][2QQ9l ). 



dMuno et al.||2006[ [Ritchie 



We note that only the two youngest SGRs have a star 
cluster in their vicinity. The lack of a star cluster in 
the vicinity of the older SGRs (despite ages of 4 to 10 
kyr) suggests that it is not essential that SGRs should be 
associated with star clusters. Furthermore, the inferred 
progenitor masses of SGR 1900+14 does not compel us to 
believe that SGRs arise from massive stars. We conclude 
that binarity likely has a bigger role in forming SGRs. 

5.2. Braking Index 

If the association of the SGRs with the star clusters is 
taken for granted, we can constrain the braking index of 
the magnetars. The braking index n is calculated from 
the following implicit equation: 

n = l + ^(l-(Po/P)("-i)). 

Here, T is the kinematic age of the magnetar (time taken 
to move from cluster to present position) and Pq is the 
spin period at birth. 

The instantaneous P is known to vary by a factor of 
three to four corresponding to la rge variations of br aking 
torque on the magnetar (Woods et al.| 
use the X-ray timing measurements ii 
[ et al.l (1998); Mereghetti et al. (2005b); Woods 



20021 12007[). We 
rom iKouveliotou 



^^;|Ma£sden et al..C1999J;. Woods et al. (2002t 
Mereghetti et al.[ ( [2UU6| r [iNakagawa et al.[ ! [2UU9p 



et al 



to cal- 
for SGR 1806-20 

and 17 x 10"^^ s s-^ for SGR 1900+14 from 1996 to 2006. 



culate an average P of 49 x 10 ss 



+0.65 
-0.24 



Assuming Pq/P <C 1, we estimate n to be 1.76 
for SGR 1806-20 and 1.16t^;^^ for SGR 1900+14. This 
is significantly smaller than the canonical value of n = 3 
for the magnetic dipole spindown mechanism for pulsars. 
Low braking indices have bee n discussed in the conte xt of 



twisted magnetospheres (eg. jThompson et al.|2Q02 ) and 
particle wind spindown (e.g. Long et al.|2U12p . However, 
the large variations in P over tens of years implies that 
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these measurements cannot be taken at face value. 
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